Optimizing the vertebrate vestibular semicircular canal: could we balance any better? 
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The fluid-filled semicircular canals (SCCs) of the vestibular system are used by all vertebrates 
to sense angular rotation. Despite masses spanning seven decades, all mammalian SCCs are nearly 
the same size. We propose that the SCC represents a sensory organ that evolution has 'optimally 
designed'. Four geometric parameters are used to characterize the SCC, and 'building materials' of 
given physical properties are assumed. Identifying physical and physiological constraints on SCC 
operation, we find that the most sensitive SCC has dimensions consistent with available data. 



Simple physical principles and scaling arguments have 
been remarkably effective in understanding organismic 
and evolutionary properties of the animal kingdom 0. 
Sensory organs, which transduce physical stimuli into 
neural signals, are particularly well-suited for physical 
analysis. The biophysics of the cochka, for example, has 
been the focus of much recent study |2| • 

Here we focus on a different inner-ear organ: the fluid- 
filled semicircular canals (SCCs) that each of the roughly 
45,000 vertebrate species employs to mechanically sense 
rotation 0, . Rotation sensation is obviously impera- 
tive for balance, but also plays an equally important role 
in vision. The neural output of the SCCs feeds directly 
to the oculomotor system, causing a reflexive motion of 
the eyes that compensates for head motion. This allows 
you to read this article even while shaking your head, 
a task which would be much more difficult if the article 
itself were shaken. 

Previously, we studied 'top-shelf vertigo', a mechan- 
ical disorder of the human SCC 0. Here, we address 
a rather astonishing feature of the SCCs: the SCCs of 
every mammal, from mice to whales, are essentially the 
same size (Fig. 0) 0,0. Mammals span seven decades in 
mass and almost three in length, yet SCC dimensions are 
restricted to less than one. In fact, the SCC of the human 
fetus reaches its full adult size by the fourteenth week of 
pregnancy 0. Fish, reptile, amphibian, and bird SCCs 
are of similar size [6j, the one apparent exception being 
sharks (discussed below) • In exploring reasons for the 
constancy of cell size, Vogel writes, "...anything biological 
that doesn't vary in size ought to strike us as notewor- 
thy" 0. Similarly, SCC non-scaling suggests that these 
particular dimensions may be special to, and perhaps op- 
timal for, their function. 

The idea that biological structures are somehow 'opti- 
mal' has recently found rather dramatic support. Micro- 
cavities in the brittle star skeleton act as perfect lenses 
0, and sea sponges develop single-mode optical fibers 
that rival current technology |lfj j. The human visual 
system operates at the single-photon level [ill, and the 
auditory system is limited by thermal noise |l2j. In fact, 
evolution can be viewed as a gradient search seeking to 
optimize 'fitness'; however, the utility of this approach 
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FIG. 1: Measured values of 7? and d for mammals 0, 0. 
Despite seven decades of mass variation, variations in R and 
d are limited to about one decade. Solid (dashed) horizon- 
tal lines reflect 'optimal' ('tolerant') SCC dimensions 1121 
1151 . Rough measurements of c and t for mammals llH reveal 
t ~ 0.15 — 0.44 mm and c ~ 0.17 — 0.5 mm (c.f. optimal 
c = 0.12, t = 0.12.) 



is often limited by the difficulty of defining 'fitness'. By 
contrast, the semicircular canals have a well-defined pur- 
pose, and reasonable measures of their quality can be 
proposed. 

Herein, we suggest that the 'universal' size of the SCCs 
can be understood in terms of a system which has been 
'optimally designed' by evolution. Others have examined 
the relation between SCC sensitivity and geometry from 
a purely fluid standpoint 0, , but did not address the 
elastic membrane and sensory hair cells that complete 
the transduction process, and which depend on SCC di- 
mensions. Our approach is as follows: we start with SCC 
'blueprints' without dimensions (Fig. |5J), and assume ba- 
sic building materials (solid walls, fluid endolymph, elas- 
tic material, and hair cells) to be given. Identifying phys- 
ically and physiologically reasonable constraints, we find 
a unique set of well-constrained and robust 'optimal' di- 
mensions that maximize SCC sensitivity. Moreover, they 
are consistent with measured data (Fig. 

We begin with a description of basic SCC structure 
(Fig. |2l . Each car contains three mutually orthogonal 
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SCCs to span the three rotation axes. Each SCC is a 
hollow torus of major radius R that consists of a narrow 
duct of radius d and a bulbous ampulla of radius c, all 
filled with water-like endolymph of density p and viscos- 
ity p. The cupula, a mucus membrane of thickness £ and 
radius c, completely spans the ampulla to block fluid flow 
(Fig. |2Jd) Q . Along one wall sits sensory hair cells, from 
which bundles of stereocilia (of length Xh « 50^m 0) 
project into the cupula. Within each bundle, stereociliar 
tips are linked by filaments that act as 'gating springs', 
opening or closing ion channels when stereociliar tips are 
displaced . Thus cupular deformations displace stere- 
ociliar tips, which cause the hair cells to fire. 

Basic SCC function is shown in Fig. An impulsive 
angular acceleration a(t) = flo6(t—t a ), viewed in the rest 
frame of the canal, gives rise to ballistic endolymph mo- 
tion with initial angular velocity f2o- Viscous resistance 
from the walls cuts off this inertial motion after a time 
Tf (~ 5 ms in humans), during which a volume V c oc Qq 
of fluid passes through the canal. The cupula distorts 
to displace the same volume (Fig. Et(i)), which deflects 
embedded stereocilia and triggers a corresponding neural 
signal. The elastic restoring force of the distorted cupula 
drives fluid back through the narrow duct. Viscous resis- 
tance from duct walls slows cupular relaxation to a time 
scale t s (4-7 s in humans and primates 0,0] ), so that V c 
encodes angular velocity VL for t < r s . Under sustained 
rotation, however, the cupula relaxes (Fig. OJi(ii)) and 
'forgets' the constant rotation rate fig. Upon stopping, 
the sudden deceleration leads to a cupular displacement 
in the opposite direction (Fig. Et(iii)), causing dizziness 
(the false sensation of rotation) . The shorter r s , the more 
quickly an organism becomes dizzy. 

Simple arguments give scaling relations for SCC pro- 
cesses. Fluid 'knows' that walls are accelerating only 
after vorticity created at the walls diffuses (with diffusiv- 
ity v = p/p) to the duct center, giving an inertial time 
Tf ~ d 2 /v. During t/, fluid moves ballistically with ve- 
locity floR, so that a fluid volume V c ~ (ird 2 )^loRTf trav- 
els through the duct. Once inertia is cut off, the cupula 
relaxes quasistatically, as the elastic restoring pressure 
AP = KV C drives a Poiseuille flow V c ~ APd 4 / p(i d R, 
assuming viscous resistance to be dominated in the nar- 
row duct. An ODE for V c results, whose solutions decay 
exponentially on a time scale t s ~ p(3 d R/ Kd 4 . 

Prefactors were obtained by solving for the time- 
dependent flow due to an impulsive acceleration, giving 



4ird 4 Ra 



^0, Tf 



8p(3 d R 
Knd 4 ' 



(1) 



were Ao « 2.4 is the first zero of the Bessel function Jo, 
and a « 1.3 is a geometric factor [lij . 

The final physical ingredient involves elastic deforma- 
tions of the cupula, which we treat as a clamped plate 
of modulus E = E'(l — a 2 ) and Poisson ratio a [l!j. A 
plate of bending rigidity Dq = E't 3 /12 obeys DqW 4 w — 




- (b) 





\Cupula\ 




V, 




2c 


Stereocilia— -f 







1 




trf 



FIG. 2: (a) The semicircular canal is a torus of major radius 
R, with a long narrow duct of radius d spanning an angle 
(3d- (b) The ampulla contains an elastic cupula of thickness 
t and radius c, which distorts under SCC rotation, causing 
embedded hair cell stereocilia to deflect a distance £ and trig- 
ger a neural signal. Measured values (in mm) for humans are 
{R h ,c h ,d h ,t h } = {2.8,0.44,0.19,0.31} 0. (c) SCC opera- 
tion: (i) Angular velocity is impulsively accelerated to Qq at 
t a and the cupula distorts by V c oc f2o. (ii) Under sustained 
rotation, the cupula relaxes and 'forgets' Qq. (iii) Upon decel- 
erating to rest at td, fluid inertia causes a (negative) cupular 
displacement, eliciting dizziness. 



AP, giving a deformed profile w = u> max [l — (r/c) 2 ] 2 , 
where w max = APc 4 /64Z?o H3- ^ n integration gives 
the cupular stiffness K relating AP C to V c , 
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whose value K = 4.6 GPa/m 3 is derived from measured 
values of r s w 4 s and QJ. As a check, this implies 
£'s;2 dyne /cm , which is consistent with measurements 
for mucus |2l| . 

Having described the mechanical response, we now 
connect cupular displacement w to stereociliar tip dis- 
placement C (and thus to a neural firing rate /). Hair 
cells have a high resting discharge rate f r and are bi- 
directionally sensitive, so that deflections in either direc- 
tion can be sensed quickly The gating springs that 
link stereociliar tips cause ion channels to open or close 
as the tips arc deflected, increasing or decreasing / lin- 
early with ( for small deflections [ljjj. Since the tips are 
embedded in the cupula, their deflection £ varies with w 
via 
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Measurements on 5-10 pxa utricular and cochlear hair 
cells indicate a threshold displacement £ m i n ~ O (nm) for 
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neural activity, and gatin g sp ring models predict Cmin to 
vary with hair cell length |l5j . Because SCC hair cells are 
about 10 times longer, one expects Cmin ~ O(10nm). A 
similar value for Cmin can be obtained from other exper- 
iments under the assumption that (order of magnitude) 
hair cell properties are conserved across species. Phys- 
iological measurements indicate a lower limit for sensa- 
tion in humans, giving a sensation threshold O m i n ~ 2°/s 
[2^ |. In squirrel monkeys, the relationship between SI (or 
C) and firing rate was measured to be linear: A/ « afi, 
where a = 0.24/° The firing rate has natural lim- 

its: / can not be negative, and fires at a maximum (in- 
jured) rate (~400 Hz). Saturating deflections Csat occur 
for A/ ~ f r w 100 Hz, implying a 'saturating' rotation 
^sat = fr/a w 420°/s. Using (|lj and ©, we find a 
dynamic range of stereocilia deflections 



20 nm, Csat ~ 4 /im, 
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and note that Cmin thus obtained is consistent with pre- 
dictions of the gating spring model. 

The sensitivity S of an arbitrary SCC follows from 
and J2J, with minimum detectable rotation Sl m in, 
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We have now characterized the mechanical-neural sig- 
nal transduction process, and can examine the effect of 
varying SCC dimensions. We scale SCC variables by 
their human values, denote scaled variables with hats, 
and pose the central question: Given 'blueprints' (Fig. 
12} and basic building materials (/1, p, E' , Xh, Cmin and 
Csat), what dimensions {R,c,d,t} maximize SCC sensi- 
tivity S = d 4 R/c 4 ? 

Certain constraints, however, must be obeyed for a 
physically- and physiologically-viable SCC. First, ther- 
mal fluctuations of the cupula must be small enough to 
escape detection, or else balance and vision would be im- 
paired. The probability of a cupular displacement fluctu- 
ation V c is given by P(V C ) ~ cxp(~KV 2 /2k B T). We use 
(0 and (0 to express P(V C ) in terms of C, giving P(C) = 
2 v / ^exp(-crC 2 ), where a = nE'c 2 t 3 /18k B TX%. The 
probability of sensing a thermal fluctuation is given by 
P(C > Cmin) = 1 - crf(cr 1 / 2 C min ). Requiring P(C > Cmin) 

gives 
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A = 0.15. 
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Note that A is rather insensitive to e: A(10 6 ) = 0.16 
and A(10~ 4 ) = 0.14. 

Second, the cupula should encode angular velocity for 
as long as possible. The smaller r s , the sooner angular 
velocity is 'forgotten' and dizziness ensues. Requiring r s 
to be larger than some minimum time t™ gives 
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FIG. 3: A 2D slice of 4D parameter space showing the con- 
strained space available for viable SCC design. Constraints 
1611 It are labeled A-F. Also plotted is the sensitivity gradient 
and, for reference, dimensions of the human SCC. The 'opti- 
mal' SCC is found at the vertex of constraints A, B, D, and 
F. 



and we take r" 1 = 1 s as a basic estimate. 

Third, the maximum cupular deflection w max for phys- 
iological rotations should be small (compared to t), to 
avoid damaging the cupula and for a linear response. The 
saturating rotation il sat is one example, giving 
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Rotations il max which saturate the neural signal, but are 
still physiologically relevant, also impose a constraint. As 
a basic estimate, we take f2 ma x = 27r rad/s, which gives 
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Finally, the ampulla must fit in the canal, 
R/c > 2c h /R h =E= .32, 
and the cupula should be plate-like for uniform C, 
c/t > tfc/ch = F = 0.70. 
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Eqs. (|5I11(I represent a nonlinear multidimensional op- 
timization, but an equivalent linear optimization is ob- 
tained by taking the logarithm of each equation. As in 
linear programming, extrema are found at vertices of the 
constraint equations (or -co, since logs can be negative). 
Here, the most sensitive SCC geometry is found at ver- 
tices A,B,D, and F, giving dimensions 
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0.38, (13) 
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and a sensitivity S 1 , 

§ = DA~ 3/5 F- 8/5 pa 24, 



(14) 



that is about 24 times greater than the human SCC. As 
seen in Fig. ^ these dimensions compare quite favorably 
with available data — particularly since such an optimum 
need not exist in the first place, and even if it did, its 
dimensions could have differed from those found in nature 
by many orders of magnitude. 

Two parameters (O roa x and r™) were chosen some- 
what arbitrarily, and it is natural to question how these 
choices affect the above results. Only R and d depend 
on these choices, and weakly at that: R ~ (rf/fimax) 1 ' 2 
and d ~ {r™Sl ms . x )~ 1 / 9 ' . For R to change by an order 
of magnitude, r™/Q max must be off by two. This inscn- 
sitivity to the crude choices made above reinforces the 
robustness of this result, and implies that existing inter- 
species variations in t™ and f2 ma x should lead to only 
slight differences in optimal SCC dimensions. 

While the above analysis assumes perfect 'machine 
tools,' nature's machinery has a rather larger tolerance: 
variations in human SCC dimension are of order 10% 
|23j . To allow for imperfect 'machining', we find the most 
sensitive 'tolerant' SCC such that 10% variations in any 
dimension satisfy constraints A — F, giving 



R = 1.9, c = 0.33, d = 0.46, t = 0.38, 



(15) 



with a relative sensitivity S ~ 7. SCC sensitivity varia- 
tions due to size variations range from about 0.4 to 2.5 
times that of the 'average' canal. 

No constraint on measurement time Tf was introduced, 
although it plays an important role in vision: to be 
'watched', an image must be kept within one retinal fovea 
- about 1° in humans In the time Tf following an Qq 
impulse, the eye/fovea lags the image by A6* ~ QoTf, 
so that the maximum rotation allowing a fixed gaze is 
flu ~ 1°/ T / ~ 200°/s. Larger creatures typically undergo 
slower rotations, allowing larger t/ (and thus d). Further- 
more, rapidly-flying birds have smaller d than similarly- 
sized earthbound ones 0] , which may allow them to hold 
objects in view more easily during rapid maneuvers. 

Although we have explicitly discussed mammals, the 
same SCC dimensions are found in all vertebrates, with 
a few notable exceptions. Head size is an obvious and 
unavoidable constraint for exceedingly small fish larvae, 
whose SCCs start small and grow with the fish 0. As 
a group, sharks-even small ones-have abnormally large 
SCCs (R ~ 40 mm) Q. Since sharks broke rather early 
from other vertebrates in evolution, a possible explana- 
tion could involve different 'building materials'. 

In summary, we have argued that evolution has con- 
verged on an 'optimal' design for a maximally sensitive 
rotation detector. The optimal SCC is well-constrained, 
with little room for variations, and falls within a factor 
of three of available data. The optimum itself is robust 



and depends on basic mechanics and established hair cell 
properties; the 'calibrating' assumptions made herein are 
self-consistent, but not essential for this central result. 

I am grateful to Howard Stone for introducing me to 
the vestibular system, and for many thoughtful com- 
ments and discussions. Conversations with Bill Bialck 
about foveation, with Michael Brenner about optimiza- 
tion, and with Eric Lauga about this article are gratefully 
acknowledged. 
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